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Abstract 

The  streamfunction-vorticity  equations  for  two-dimentional  cavity  flow  are  solved  by  a  new  finite 
element  method  which  uses  finite  spectral  basis  functions  as  shape  functions  for  rectangular  elements. 
Simulations  for  several  cases  with  different  Renolds  numbers  are  performed.  Good  agreement  was 
obtained  in  the  comparison  between  the  present  results  with  the  bench  mark  solutions. 
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1.  Introduction 


During  the  past  decades,  Finite  Element  Methods  (FEMs)  have  been  recognized  to  be  powerful 
tools  in  the  solution  of  numerous  flow  problems.  More  and  more  methods  based  on  FEM  were 
developed  and  yielded  many  satisfactory  results  in  various  problems. 

The  finite  spectral  method  based  on  non-periodic  Fourier  integral  has  succeeded  in  dealing  with 
spectral  methods  pointwise  [1,2]  and  has  been  successfully  applied  to  classical  schemes  such  as  NND, 
ENO  etc[3,4].  It  is  characterized  by  local  property,  non-periodicity,  orthogonal  relation,  efficiency  and 
simplicity.  This  makes  it  possible  to  apply  finite  spectral  method  to  finite  element  method  by  using 
Wang  Kemal  as  Shape  functions.  In  order  to  verify  whether  it  is  applicable,  we  take  driven  cavity  flow 
as  a  test  problem  and  compare  the  numerical  results  of  this  problem  with  the  benchmark  solutions. 


2.  Wang  Kernal 

Expanding  the  dicrete  pulse  function 


/,= 


[1  ifj=0 


0  ifj=others 

by  Fourier  series  and  truncated  off  the  Nth  term,  we  obtain 


(1) 
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where  —l<x<l,  and  Cn  takes  0.5  for  n  =  ±N  but  1  for  others.  Since  it  is  similar  to  Dirichlet 

Kemal  which  is  difined  on  the  infinite  interval,  we  call  it  Wang  Kemal.  It  is  easy  to  expand  one¬ 
dimensional  Wang  Kemal  to  two  dimensions.  Two-dimensional  finite  spectral  shape  function 
WN  ( x ,  y)  is  represented  as 


WN(x,y)  =  WN(x)*WN(y) 


Flence  WN  (x,  y )  can  be  used  as  basis  functions  in  two-dimensional  problems. 


(3) 
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3.  Computation  of  Cavity  Flow 


3.1  Problem  description 

The  computation  of  the  cavity  flow  in  the  square  domain  has  been  viewed  as  one  of  the  standard 
test  problem.  The  problem  consists  of  a  square  cavity  totally  filled  with  an  incompressible  viscous 
fluid  and  a  top  wall  moving  with  constant  velocity.  The  problem  description  is  shown  in  Fig.  1. 
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Fig.  1.  Driven  cavity  flow  problem  description 


3.2  Governing  equations 

For  two-dimensional  incompressible  laminar  flows,  the  Navier-Stokes  equations  can  be  written 
in  streamfucntion-vorticity  formulation.  The  streamfimction  and  vorticity  equations  can  be  written  in 
dimensionless  form  as  [5] 


d2y/ 

8Xf  dx[ 


=  -co 


(4) 


and 
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where  (//  is  the  streamfunction,  X;.  (i  =  l,  2)  are  dimensionless  Cartesian  co-ordinates, 


y  = 


8y/ 

ay 


(6) 


and 


U2  = 


8y/ 

ay 


are  the  velocity  components  in  the  Xx  -  and  X2  -  directions  respectively, 


ay  ay 

(O  =  — - - - 

ax,  ax2 


(7) 


(8) 
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is  the  vorticity,  t  is  the  time  and  Re  is  the  Reynolds  number. 

The  boundary  conditions  of  driven  cavity  flow  are  given  as  follow,  the  walls  are  no-slip 
boundaries,  the  value  of  the  streamfunction  is  known 


¥  =  0  (9) 

since  the  walls  are  also  streamlines.  In  the  procurement  of  the  wall  vorticity,  utilize  the  tangential 
component  of  velocity  and  the  streamfunction’s  value  near  walls[6], 

3.3  Finite  element  formulation 

In  the  elment,  the  unknown  variables  can  be  approximated  by  means  of  the  standard  expansions 


n 


¥  ~  Z  Nixffi 

i= 1 

(10) 

n 

NfiJi 

i= 1 

(11) 

where  Nt  are  two-dimensional  finite  spectral  basis  function,  I//  and  coi  are 

nodal  values  of  y/  and  co 

respectively  and  n  is  the  number  of  nodes  in  the  element.  Following  the  Bubnov-Galerkin  method, 
equations  (4)  and  (5)  can  be  written  in  matrix  form  as 

Aiyt-Bfa,  =  Cf 

(12) 

(13) 

In  the  above  equations 

4’=l Lj»,x,NJ.x+Nl,xN.Xi)dn 

(14) 

B?  =  \LN'N>dn 

(15) 

cr'=\i„NMrdn 

(16) 

4‘’=Lw« 

(17) 

C  =JL N, dn 

(18) 

Ff' N,,*,wn 

(19) 

(20) 

q  and  qm  are  natural  boundary  conditions  of  (//  and  CO  respectively. 

3.4  Results  and  discussion 

Using  finite  element  method  with  finite  spectral  basis  functions,  the  streamfunction-vorticity 
equations  for  two-dimentional  cavity  flow  are  uncoupled  and  solved  in  sequence  on  uniformly  spaced 
grids  of  elements  consisting  of  21X21  mesh  points.  The  discrete  elements  are  4-node  rectangular 
elements.  Results  for  Re  =  100, 400, 1000  are  obtained.  The  figures  of  contours  of  streamfunction 
and  contours  of  vorticity  for  Re  =  400  shown  in  Fig.  2  are  similar  to  previous  work  [5]. 
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Fig.  2.  Contours  of  streamfunction  and  contours  of  vorticity  for  Re=400 

In  Fig.3  and  Fig.4  we  compare  the  centerline  velocity  profiles  with  benchmark  solutions  [7].  the 
results  are  not  in  well  agreement  When  Re=1000.  Consequently  we  increase  mesh  points  to  31X31 
and  acquire  better  results  shown  in  Fig.  5. 


Crvclocity  along  vertical  centerline  {/^-velocity  along  horizontal  centerline 

Fig.  3.  Velocity  profiles  along  the  centreline  with  21 X21  mesh  at  Re=100 


u,  X, 


C'l-vclocity  along  vertical  centerline  FVvcIocity  along  horizontal  centerline 

Fig.  4.  Velocity  profiles  along  the  centreline  with  21 X21  mesh  at  Re=1000 
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(Ti-velocity  along  vertical  centerline  CVvelocity  along  horizontal  centerline 

Fig.  5.  velocity  profiles  along  the  centreline  with  31 X31  mesh  at  Re=1000 


In  Fig.  6.  we  compare  the  centerline  velocity  profiles  with  previous  finite  element  solutions[8] 
and  find  out  that  They  are  in  well  agreement . 


Fig.  6.  comparison  with  previous  finite  element  solution[8] 


4.  Conclusions 

In  this  paper,  we  use  finite  spectral  basis  functions  as  interpolation  functions  for  rectangular 
elements  and  solve  the  streamfiinction-vorticity  equations  for  two-dimentional  cavity  flow.  Results  are 
in  well  agreement  with  previous  solutions.  These  prove  that  the  spectral  basis  functions  are  applicable 
and  effective  and  this  method  can  be  used  in  solving  other  problems. 
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